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Recently observed increase of direct CP asymmetry in charm meson nonleptonic decays is difficult to explain 
within the SM. If this effect is induced by new physics, this might be investigated in other charm processes. We 
propose to investigate new CP violating effects in rare decays D — > P£ + £~ , which arise due to the interference 
of resonant part of the long distance contribution and the new physics affected short distance contribution. Per- 
forming a model independent analysis, we identify as appropriate observables the differential direct CP asym- 
metry and partial decay width CP asymmetry. We find that in the most promising decays D + — > ir + £ + £~ and 
Df — ¥ K + £ + £~ the "peak-symmetric" and "peak-antisymmetric" CP asymmetries are strong phase dependent 
and can be of the order 1 % and 10 %, respectively. 
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I. INTRODUCTION 



In last two decades chances to observe new physics in charm processes were considered to be very small. In the case of 
flavor changing neutral current processes the Glashow-Iliopoulos-Maiani (GIM) mechanism plays a significant role, leading to 
cancellations of contributions of s and d quarks, while intermediate b quark contribution is suppressed by V u b element of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix. The lack of appropriate theoretical tools to handle long distance dynamics in these 
processes is even more pronounced than in the case of B mesons due to abundance of charmless resonances with the masses 
close to the masses of charm mesons. However, this has changed at the end of last year when LHCb experiment reported a non- 
vanishing direct CP asymmetry in D° — > K + K~ and D° — > tt + ii~ [1] also confirmed by the CDF experiment [2]. Many papers 
investigated whether this result can be accommodated within the standard model (SM) or is it new physics (NP) that causes such 
an effect. The measured difference between the CP asymmetry in D a — > K + K~ and D° — > it + tt~ is a factor 5—10 larger 
than expected in the SM and eventually can be a result of nonperturbative QCD dynamics as pointed out in refs. [3-10]. Model 
independent studies [11, 12] indicated that among operators describing NP effect, the most likely candidate is the effective 
AC = 1 chromomagnetic dipole operator. In order to distinguish between SM or NP scenarios as explanation of the observed 
phenomena it is crucial to investigate experimentally and theoretically all possible processes in which the same operator might 
contribute. Recently the effects of the same kind of new physics have been explored in radiative [13] and inclusive charm decays 
with a lepton pair in the final state [14]. In [13] it was found that NP induces an enhancement of the matrix elements of the 
electromagnetic dipole operators leading to CP asymmetries of the order of few percent. 

In addition to radiative weak decays, charm meson decays to a light meson and leptonic pair might serve as a testing ground 
for CP violating new physics contributions. As in other weak decays of charm mesons the long distance dynamics dominates 
the decay widths of D — > P£ + £~ [15-17] and it requires special task to find the appropriate variables containing mainly short- 
distance contributions. In this study we investigate partial decay width CP asymmetry in the case of D — > P£ + £~ decay. The 
short distance dynamics is described by effective operators 07, 0g, and 0io of which the electromagnetic dipole operator O-j 
carries a CP odd phase of beyond the SM origin, developed due to mixing under QCD renormalization with the chromomagnetic 
operator. In this paper we investigate impact of this mixing on the D — > P£ + £~ decay dynamics. The paper is organized as 
follows: Section II contains the description of the short distance contributions, Sec. Ill is devoted to the long distance dynamics. 
In Sec. IV we present the partial width asymmetry. We summarize our findings in Sec. V. 
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II. CP-ODD EFFECTIVE HAMILTONIAN 

The genuine short distance (SD) dynamics of c —> u£ + £~ decay on scale ~ m c is defined by the effective Hamiltonian 

V 2 i=7,9,10 
6/71 

^7- t^^PbcF"", 



(4tt) 5 

e 2 
(4tt) 5 

e 2 



9 = 7 j^(urP L c)(h^). 



SD amplitude has contributions from electromagnetic penguin Or and four-fermion operators with vector and axial lepton 
current, Og and 1O , respectively, and while their contribution to the decay width is negligible in the resonance-dominated 
regions, possible imaginary parts of Wilson coefficients may generate direct CP violation via interference with the CP-even 
long distance amplitude. Here we will assume that CV(m c ) has a large NP phase. It can be naturally induced by QCD mixing 
between CV and C$ and an assumption that direct CP in singly Cabibbo suppressed decays is explained by new physics in C%. 
We consider the range proposed in [13], 

|ImCV(m c )| = (0.2 - 0.8) x 10" 2 . (2) 

Unlike in the B meson decays, contributions of the third generation quarks in the flavor changing neutral current loops are 
not favored over the first and second generation quarks in charm decays. Effects of virtual light quarks cannot be expanded in 
series of local operators as in (1) since they are genuine long-distance (LD) effects and we will address them in the next section. 
Keeping only the CP-violating part of the SD amplitude of D — > nl + t~, where I = e, //, we have 

„cpv iV2G F a m c 2 - 

-4sd = C 7 {m c ) ■ fT{q){ip£), (3) 

where p is momentum of the D meson. The form factors for D — > it transition via vector current and electromagnetic dipole 
operators are defined as customary 



(n(p')\u lfl c\D(p)) 



(P + P)m~ ,2 ^ 



2 2 

m n — mi 



Ji(g a )+ ° 2 "N ^ (g 2 ). 



q 



(7T(p')\ua^c\D(p)) = -i (pX - p v p^) 2fr( f 2) , (4) 

with q 2 = (p — p') 2 . In the heavy quark limit of large m c the matrix elements of tensor and vector operators are related in the 
rest-frame of the heavy meson [ 1 8] 

(tt( P ') I ua i0 c | D(0)) = t (ir(p') \ u lt c \ D(0)) . (5) 
The form factor Jt is then connected to F\ and F [19] 

2fr^)^- = F l{q 2 ) ( 1 + - F Qiq 2 )^t . (6) 

m D +rrin \ q J q 

Close to the kinematical point of pion at rest, q 2 as m D , the form factors scale as F\/Fq ~ mo [18, 20] and the above 
relation reduces to fx ~ F\. Form factors Fi (q 2 ) for D —> it transitions have been computed in unquenched lattice QCD 
simulation with staggered light quarks in ref. [21], There the shapes of the form factors are presented in the Becirevic-Kaidalov 
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parameterization [22], 



F ^ 2 ) = / 1 _^w° 1 ) _ o _ g _ V (7) 



Fi(0) 

_ 1 q 2 
b TnT 



Fi(0)= 0.64(3), (8) 
a = 0.44(4) , 
6= 1.41(6). 

We employ the form factor /t(<7 2 ) as given by heavy quark symmetry relation (6) used with shapes (7), (8) for Fi(q 2 ) and 



III. RESONANT AMPLITUDE 



The <p resonant peak in the lepton invariant mass distribution of D —> P£ + £~ is very narrow (T^/rriA ~4x 10~ 3 ) and well 
separated from other vector resonances. Therefore, close to the peak the long distance amplitude is, to a good approximation, 
well dominated by a Breit-Wigner amplitude [15-17] 

A* D [D -> ncf> -> n£-£+] = ^V^vJ^a^* 2 ^ *{kJ)fv{k + ) ■ (9) 

Here we use a = 1/137 in the leading order in electromagnetic interaction. Finite width of the resonance generates a q 2 - 
dependent strong phase that varies across the peak. We have also introduced the strong phase on peak, 5$, and the normalization, 
a^, that are both assumed to be independent of q 2 . Parameter is real and can be fixed from measured branching fractions of 
D — > tt<P and <j> — > £ + £~ decays [17]. For definiteness we will focus on the I = /x decay modes. From the Particle Data Group 
compilation we read [23] 

Br(L>+ ->■ 07T+) = (2.65 ± 0.09) x 10~ 3 , (10) 
Br(</> -> fj, + fj,~) = (0.287 ± 0.019) x 10~ 3 , 

and when we take into account the small width of <f» 

Br(D + 7r + 0(^ m + aO) « Br(L>+ -> 0tt+) x Br(0 ->■ /J+/0 , (11) 

we find from eq. (9) 

a = 1.23 ±0.05. (12) 



IV. DIRECT CP ASYMMETRY 



The direct CP violation in the resonant region is driven by the interference between the CP-odd imaginary part of the SD 
amplitude and the LD amplitude. The pair of CP-conjugated amplitudes read 

A{D+ ^^l+n^A^ + A^ , (13) 
A(D-^n-£+r) = Ai Ti + As™ . 

In principle the short-distance amplitude contains a strong phase that can be rotated away because the overall phase of the total 
amplitude is irrelevant. The CP-odd part of the LD amplitude is proportional to the imaginary part of the relevant CKM factor 



V* S V US = A + iA 2 v \ 5 , 



(14) 
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and can be safely neglected and accordingly we have A LD = A LD . Then the differential direct CP violation reads 

\A 2 - W 



acpWq 2 ) 



|-4| 2 + |-4| 2 

-3 f T (q 2 ) 



(15) 



ImC 7 



27r 2 m D + m v Re[V* s V us ] 



COS <L 



m^T 4 



sin 5 d 



When considering numerics in what follows we will set ImCV to the largest possible value of 0.8 x 10~ 2 for the sake of 
illustration. Relative importance of the cos 5^, and sin 6^ for representative choices of S^, is shown on the upper plot in fig. 1 . The 
linearly rising behaviour of the sin ^-driven term of the asymmetry is compensated by a rapid drop of the resonant amplitude (9) 
that severely diminishes number of experimental events as we move several away from mu — rriip- Both effects are included 
in the effective experimental sensitivity that also takes into account the rate of events in the considered kinematical region and is 
shown on the bottom plot of fig. 1. There we plot acp{mu), weighted by the differential branching ratio, a combined quantity 
that scales as ~ A^lmAg^ . These sensitivity curves expose entirely different behaviour than acp(m«). If the phase 5<p 
is close to or ir one finds the best sensitivity close to the peak. On the contrary, for 5$ ~ ±7r/2, the CP asymmetry is an 
odd-function with respect to the resonant peak position and is maximal when we are slightly off the peak. Therefore, experiment 
collecting events in a symmetric bin around m« = would be unable to observe CP asymmetry for maximal phase S ~ ±7r/2. 



A. Partial-width CP asymmetries 

In order to keep the experimental search as general as possible one should use appropriate search strategies to address the 
two limiting possibilities, i.e. <5^ = 0,tt and 8^ = ±ir/2. First, let us define a CP asymmetry of a partial width in the range 

mi < mu < mi, 

r(mi < m u < m 2 ) - f (mi < m u < m 2 ) 

Acp{mi,m 2 ) = — N r, (16) 

1 (mi < mu < m 2) + 1 \jni < mu < m 2 ) 

where T and T denote partial decay widths of D + and D~ decays, respectively, to 7r ± /i + /i _ . Acp is related to the differential 
asymmetry acp ( Vy) as 

Acp{m 1 ,m 2 ) = -a , (17) 



where 

i /-s m ax(<? 2 ) . 

g2_ m 2 )2 + TO 2 r 2 / | if 



(18) 



involves the resonant shape and the integral of the lepton trace over the Dalitz variable s = (p' + fc_) 2 whose kinematical limits 
read 



(q 2 ^l- 4 -^T\ 1/2 (q 2 ,ml,ml)y 



*max/minV« ) . 2 ^ 2 ' yLy > 

X(x, y, z) — (x + y + z) 2 — 4(xy + yz + zx) . 

The D + — > 7r + e + e~ decay mode been searched for by the CLEO experiment [24] where signal in a bin around the <f> 
resonance was observed. The following partial branching ratio was reported 

Br(£>+ -> 7r+e+e-)| mcc _ m ^|< 20 Mo v - (1.7 ± 1.4 ± 0.1) x lO" 6 , (20) 

in a bin up covering the region ~ 5 T$ to the left and right from the nominal position of the (j> resonance. We define the 
asymmetry on same bin for the 7r + ^ + /i~ final state as 

C^p =A C p(m - 20 MeV,m + 20 MeV) . (21) 
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FIG. 1. Top: CP asymmetry acp(mu ) around the <f> resonance (dashed vertical line) for representative values of strong phase 5^ — 0, n/2, n. 
Bottom: (dBr/dmii) acp(mu), the measure of sensitivity to direct CP-violation. Dashed vertical lines at mu — ± F^, denote the width 
of the resonance. 
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The asymmetry Cq P is most sensitive to the cos 5^ term in Eq. (15) and is therefore optimized for cases when 5$ ~ or 6$ ~ 7T. 
Its sensitivity would decrease if we approached 8$ ~ ±7r/2, since the acpijnti) would be asymmetric in (m« — m^) in this 
case. For that very region of 6a, we find the following observable with good sensitivity to direct CP violation 

S£ p = Acpimj, - 40 MeV, - 20 MeV) - A CP (m + 20 MeV, m + 40 MeV) (22) 

The bins where the partial width CP asymmetries Cq P and Sq P are defined are shown in fig. 2 together with acp(mu). 




0.98 1.00 1.02 1.04 1.06 0.98 1.00 1.02 1.04 1.06 

m n [GeV] mu [GeV] 



FIG. 2. Left: Asymmetry acpirau) weighted by dBx/dmu m me case when dominated by cos 5^ term. The shaded region denotes the 
defining bin for asymmetry C^ p . Right: acp(jntl) when dominated by sinj^. Shown are also the two bins where the asymmetry Sq P is 
defined as the difference of Acp in the two bins. 



B. Case study for C*^ p and S$ P 

The asymmetry Sq P can be an order of magnitude bigger than Cq P (see fig. 3, left). However, when we rescale the asymme- 
tries by the branching ratios in the bins where these asymmetries defined, namely by 7.1 x 10~ 7 for Cq P and 6.7 x 10~ 8 for 
Sq P , we find evenly distributed sensitivity to direct CP violation over entire range of 5$. Also in the transient regions between 
the regimes where either cos 6$ or sin 5^ terms dominate the sensitivity does not decrease significantly. Numerical values of 
the central values are summarized in tab. I, whereas the errors coming dominantly from parameter (12) and the form factor 
fx (6), (8) are estimated to be of the order 20 %. 



5$ 


Cgp x 10 2 


Sgp x 10 2 


Br(C-bin) C^ P x 10 7 


Br(S-bin) Sf, P x 10 7 


0,7T 

±tt/2 


T0.26 
±0.005 


±0.014 
T6.8 


=F0.018 

±4 x 10~ 4 


±4 x 10" 5 
=F0.017 



TABLE I. Values of D — > 7r + /i + /i CP asymmetries Cq P and Sq P for representative values of 5$. Last two columns show effective 
sensitivity. 



C. Comment on D s -> <f>K + K + £ + £~ 

Same type of asymmetries can be defined for the decay mode of D s meson via the (f> resonance to final state K + £ + £~ . The 
resonant amplitude is described by an analogous expression to (9) and is parameterized by real aL and 5'^. The branching ratio 

Br(D s -> cj)K + ) = (1.8 ±0.4) x 10~ 4 , (23) 
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FIG. 3. Partial width asymmetries of D — ¥ n + £ + £~ decay. Left: asymmetries Cq P and Sq P for ImCV = 0.8 x 10~ 2 and their dependence 
on 5$. Right: asymmetries rescaled by the branching ratios in the corresponding bins, thus representing effective sensitivity to direct CP 
violation. 



is an order of magnitude smaller than the corresponding Br(D + — > cf>ir + ). By employing the narrow width approximation the 
value we find a'^ = 0.49 with ~ 10 % error. On the other hand, the short distance amplitude remains of same order of magnitude 
as in the D + — > ir + /i + fi~ case. We neglect the S'C/(3)-breaking corrections to the form factor and use /t(<? 2 ) as given in (6) 
adjusted by — > vtik- The asymmetries Cq P and Sq P are larger, whereas the experimental sensitivity is weaker due to smaller 




FIG. 4. Partial width asymmetries of D s — ¥ K + l + £~ decay. Left: asymmetries Cq P and S^ p for ImCV = 0.8 x 10 -2 and their dependence 
on 5's. Right: asymmetries rescaled by the branching ratios in the corresponding bins, thus representing effective sensitivity to direct CP 
violation. 

branching fractions, as shown in tab. II. 
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S', 


C^p x 10 2 


Sg'p x 10 2 


Br(C-bin) C^p x 10 7 


Br(S-bin) S^' P x 10 7 


0,7T 


T0.71 


±0.038 


T0.0035 


±2 x 10~ 5 


±tt/2 


±0.012 


Tl9 


±6 x 10~ 5 


=F0.009 



TABLE II. Values of D 3 — > K + /j, + fi CP asymmetries Cq P and Sq P for representative values of Last two columns show effective 
sensitivity. 

V. SUMMARY 

In this article we have studied CP asymmetries of rare decays D + — > n + fi + n~ and D s — > K + fi + n~ defined close to the <p 
resonance that couples to the lepton pair. These asymmetries can be generated by imaginary parts of Wilson coefficients in the 
effective Hamiltonian for c — > u£ + £~ processes. We have limited the discussion to the electromagnetic dipole coefficient C-j 
which can carry a large CP odd imaginary part, if the direct CP violation in singly Cabibbo suppressed decays D — > PP is to 
be explained by NP contribution to the chromomagnetic operator C 8 . 

We have focused on the CP asymmetry around the <fi resonant peak in spectrum of dilepton invariant mass. There approximate 
description of the resonant amplitude by means of the Breit-Wigner ansatz with two additional parameters seems valid. We have 
fixed one of these parameters from the known resonant branching fractions of D/ s ^ — »■ <f>{—¥ fi + fi~)P, while the remaining 
parameter is an unknown CP even strong phase 5$. The Breit-Wigner shape contains an additional phase that depends on the 
dilepton invariant mass. The hadronic dynamics of the short distance part of the amplitude is contained in a tensor form factor, 
fx, that can be related to F\ and F form factors using the heavy quark symmetry. The shapes of the latter two form factors 
were determined on the lattice. 

The interference term between the resonant and the short distance amplitude that drives the direct CP asymmetry depends 
decisively on the particular value of the strong phase. Namely, for large strong phase 5$, i.e., close to either ±7r/2 or — tt/2, 
the CP asymmetry would vanish if the experimental bin enclosed the 4> peak symmetrically. Conversely, the same asymmetry 
would be most sensitive to direct CP when the strong phase was either around or n. In order to cover experimentally the 
whole range of strong phase values we have devised two asymmetries that are maximally sensitive either to peak-symmetric or 
peak-antisymmetric CP violation. Taking 0.008 for the imaginary part of C\, the two asymmetries can take values of the order 
10% for 5 $ — ±7r/2 or of the order 0.1 — 1 % for 5^ = 0,n. When we multiply the asymmetries by the partial branching 
fractions in the corresponding bins, the two asymmetries provide almost even sensitivity for all values of the strong phase. 
For the D — > tt + ^ + ^ thus defined sensitivity amounts to ~ 1.5 x 10~ 9 and <~ 4.5 x 10~ 10 for D s — > K + [i + [i~. We 
conclude that measurements of partial width CP asymmetries in decays D + — > 7r + /j+/i~ and Z?+ — > K + [i + /j,~ might be 
useful in investigating whether new physics in chromomagnetic operator is responsible for direct CP violation in singly Cabibbo 
suppressed decays to two pseudoscalar mesons. 
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